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Our objective is to study resonant tunneling of an electron in the presence of inelastic scattering 
by optical phonons. Using a recently developed technique, based on exact mapping of a many- 
body problem onto a one-body problem, we compute transmission through a single site at finite 
temperatures. We also compute current through a single site at finite temperatures and an arbitrary 
strength of the potential drop over the tunneling region. Transmission vs. incident electron energy 
at finite temperatures displays additional peaks due to phonon absorption processes. Current at 
a voltage bias smaller than the phonon frequency is dominated by elastic processes. We apply 
the method to an electron tunneling through the Aharonov-Bohm ring coupled to optical phonons. 
Elastic part of electron-phonon scattering does not affect the phase of the electron. Dephasing 
occurs only through inelastic processes. 

PACS numbers: 73.23.Hk, 71.38.+i, 73.50.Bk, 03.65.Bz 



I. INTRODUCTION 

Advances in crystal-growth techniques and constantly 
shrinking semiconductor devices have motivated re- 
searchers to study electron tunneling in mesoscopic struc- 
tures. The ability to grow nearly perfect microscopic 
devices has enabled experimental confirmation of many 
theoretical predictions based on rather simple micro- 
scopic models. In particular, some basic ideas about 
electron-phonon interaction motivated an explanation of 
the seemingly unusual properties of electrons tunneling in 
the presence of interactions with phonons and other exci- 
tations , where inelastic processes affect the peak-to- 
valley current ratio. This is important in device applica- 
tions. Another line of research in mesoscopic structures is 
focused on Aharonov-Bohm oscillations in a mesoscopic 
ring (^-Q]. Special attention in this field is devoted to 
a loss of coherence, or dephasing, under the influence of 
inelastic scattering leading to suppression of hc/e oscil- 
lations. 

Most previous treatments of this problem use a Green's 
function approach often based on Keldysh formalism 
|^| O. Exact solutions were obtained only in the wide- 
band limit and for certain special cases. Calculating the 
current in the presence of inelastic scattering, some au- 
thors work in the limit of large bias through the tunneling 
region. This results in neglecting the backward current 
and the exclusion of the filled final states || . Indepent- 
boson model has been successfully used to directly solve 
the ID Schrodinger equation for arbitrary barrier struc- 
ture Transfer-matrix approach has proven very ef- 
ficient to model realistic barrier structures and compare 
calculations to experiment. Electron-phonon interaction 
is in this case treated via the Fermi golden rule |fl4) . 
Other attempts rely on linear response theory using the 
Kubo formula 0. 

This work is based on a recently developed method |ll| 
for studying inelastic electron tunneling. The method 



provides exact solutions for tunneling problems where a 
single electron tunnels in the presence of phonon degrees 
of freedom that are limited to the tunneling region. The 
main goals of this work are: a) to extend the existing 
method to finite temperatures, b) to derive an approx- 
imate formula for electric current in the presence of in- 
elastic degrees of freedom, and c) to apply developed for- 
malism to the case of a single site coupled to phonons 
as well as to an Aharonov-Bohm ring. The method we 
use allows numerically exact calculation of transmission 
at finite temperatures. We treat bands exactly in con- 
trast to some previous works ||||, where the wide-band 
limit was used. We derive an approximate formula for the 
electric current that contains an exact expression for the 
transmission matrix and provides results for current at an 
arbitrary voltage drop through the tunneling region, as 
precisely as the one-electron approximation allows. In de- 
riving the equation for current we treat the forward and 
backward current on an equal footing. We also take into 
account exclusion from the filled final states. However, 
we neglect effects of Coulomb interaction and phonon 
mediated electron-electron interaction. Our method can 
be applied to complicated tunneling structures contain- 
ing multiple- phonon degrees of freedom. In this work 
we present results for two cases a) a single site coupled 
to optical phonons and b) to an Aharonov-Bohm ring 
in a tight-binding approximation where each site is cou- 
pled to Einstein phonons. We study the effect of inelastic 
scattering on transmission and consequently the current 
through the ring. We finally comment on the effect of 
inelastic scattering on Aharonov-Bohm oscillations. 

This work is organized as follows. In Sec. II we in- 
troduce the Hamiltonian for resonant tunneling in the 
presence of inelastic degrees of freedom (optical polarons) 
that are limited to a small region of space. We describe 
the method and derive an approximate equation for cur- 
rent. In Sec. Ill we present results for transmission and 
current, at zero and finite temperature, through a single 
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site coupled to phonons. In Sec. IV we give results for 
transmission and current at zero temperature through 
the Aharonov-Bohm ring. Concluding remarks and sug- 
gestions for future work are presented in Sec. V. 

II. METHOD 

The Hamiltonian we use can be written as a sum of 
the electron part H e i, the phonon part H p h and finally 
the electron-phonon interaction H e i^ p h 

H = H e [ + H, p h + H e i-ph, 

H e l = e 3 C ) C i ~ t 3 ,k{c)c k + h.C.), 

Hph — W Y^ a rn a m, 
m 

H e l-ph = ~ Y X 3 C ]°j ( a ] + a j)- (!) 

i 

The potential, ej on site j, can describe a tunnel bar- 
rier or a voltage bias. Since we treat left and right lead 
exactly, we choose a constant potential &l(r) within the 
left (right) lead. The voltage drop is thus limited to the 
tunneling region. The hopping amplitude tj t k is set to t 
within the leads, to is the hopping amplitude from the 
lead to the tunneling region, and Xj is the (diagonal) 
coupling of an electron on site j to the phonon mode 
on the same site. Electron-phonon coupling is limited to 
the tunneling region. We consider dispersionless Einstein 
phonons with frequency to. 

The problem we are facing is to solve a scattering prob- 
lem of a single electron in the presence of inelastic degrees 
of freedom. Since a detailed explanation of the method 
was given in previous work |Iq ], we will present only a 
short overview that is necessary for the reader to under- 
stand our generalization to finite temperature and com- 
putation of the current. 




FIG. 1. Each dot represents a basis state wavefunction 
?P(j,n) in the many-body Hilbert space. The lowest row of dots 
are the sites j — — 3, . . . , 3 with diagonal energies ej . The rows 
above represent the same sites with n = 1, 2 phonon quanta on 
the site j — 0. Their diagonal energies are ej + ui and tj + 2ui, 
respectively. The bonds represent non-zero off-diagonal ma- 
trix elements in the Hamiltonian. The horizontal bonds are 
the hopping amplitudes tj t k- The vertical bonds represent the 
electron-phonon interaction. The dots can also be interpreted 
as Wannier orbitals in an equivalent 1-body tight-binding 
model. 

Consider for simplicity the case where the tunneling re- 
gion consists of a single site (j = 0) with a single phonon 
mode that couples to the electron density on site (see 
the lowest row in Fig. (M). Hopping matrix elements 
are tt,i = t for the lead nearest neighbors k, I 0, and 
tk.i = to between site and sites ±1. The wavefunction 
can be written as ip(j,n)i where the site index j represents 
the position of the electron and n represents the number 
of phonons on the site j — 0. The Schrodinger equation 
can be written for this simple problem in the compact 
form, 

- XS jfi (Vn + l^o>+i) + Vni>v,n-i)) ■ ( 2 ) 

It is already apparent from Eq. (g) that a one- 
dimensional many-body problem, consisting of an elec- 
tron and different numbers of phonon quanta, can be vi- 
sualized as an effective two-dimensional 1-body problem 
with n as the second dimension. For a better perception 
we present Eq. (^) graphically in Fig. (Q) as a tight bind- 
ing model in two dimensions. For illustrative purposes 
we have restricted the variational space to a maximum 
Nph — 2 phonon quanta. At zero temperature, an elec- 
tron incident from the left is an incoming plane wave on 
the lower left lead (there are no excited phonons on site 
0) . It has an amplitude to exit on any of the six leads, 
corresponding to elastic and inelastic backscattering and 
transmission. At finite temperatures a finite number of 
phonon quanta n are excited on site j — before scatter- 
ing with the probability P(n) = (1 - e -^) e -"^ . Thus 
an electron can enter on any of the horizontal leads corre- 
sponding to different n. We can solve the set of equations 
in Eq. (|2J) taking into account the boundary condition 
specifying that an electron can enter only through one 
lead at a time. Solutions within the leads for an electron 
approaching from the left are: 

^0<o,») = exp (ik^j) + BW cxp j) 
^0>o,m) = C (m) exp (ik% n) 3 ) 

Aj<0,m^n) = CXp HAft), (3) 

where n(m) represents the number of excited phonon 
quanta before (after) scattering, is the amplitude 
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of the incident wave, and and represent re- 

flection and transmission amplitudes. Wavevectors are 
defined by the conservation of total energy, 



+ nu — 2tcos(k^) = e l(r) + muJ — ^ C0S (^ijR)) 



(m) 



(4) 



Using the pruning technique ||15| we can remove all the 
leads that contain outgoing waves from the calculation. 
We are thus left with a system of linear equations that 
connect wavefunctions on the central site ip(o tn ) with the 
lead that carries the incoming electron. This problem can 
be solved easily by recursion for essentially any number 
of phonon quanta. The transmission matrix T L '_^ R (e, e') 
is defined as Ju| 



rp(n,m) ( /n 
1 L^R K^i fc I 



AW 



2 . j (m) 

sin k R 



sin k 



(n) 



(5) 



where n(m) represent incoming (outgoing) channels (by 
a channel we denote the lead with a specified number of 
phonons (n)), and e (e') are incoming (outgoing) elec- 
tron energies counted from the middle of the band i.e. 
e = — 2t cos(fc^)- We define the total and elastic trans- 
mission as a sum over all incoming channels n, weighted 
by the probability P(n) and a sum over all outgoing chan- 
nels. Since P(ri) depends on the temperature of the tun- 
neling region, we must assume, that the temperature of 
the tunneling region is well defined. This is achieved by 
coupling it to the external heath-bath. 



T tot (e) = ^P(n)T^(e, e ! 
T elast (e) = ^P(n)Tti(e,e'). 



(6) 
(7) 



Using the presented technique we can compute trans- 
mission at zero and finite temperatures with any desired 
accuracy. However, transmission is not in principle a 
directly measurable quantity. One has to compute elec- 
tric current in order to provide a measurable quantity. 
The total transmission as defined in Eq. (||), enters the 
equation for current only in the case of a high-voltage 
bias where neither backward current nor the exclusion of 
filled states in the right lead are taken into account ||. 

In deriving the equation for the current we start by 
observing that inelastic processes can be viewed as mul- 
tichannel electron tunneling. The current from the left 
lead, due to an electron entering the tunneling region 
through a channel n and exiting through m, is given 

(n) 

by the integral over the incoming momenta fidk L times 
the velocity of the incoming electron vl = [l/h)de/dk^£\ 
times the transmission probability T^_^ (e, e'), times the 
appropriate combination of Fermi functions from the left 
and the right lead. The total current from the left is 



expressed as a sum over all channels weighted with the 
appropriate Boltzmann factor. The net current is given 
finally by the difference between the right- and the left- 
flowing currents 

/It 
de^2P(n)Tt:^(e,e')f L (e) (1 - f R (e')) 
-it t~t 



— / de'^ 



P(m)T R ^l\e',e)f R (e')(l-f L (e)), (8) 



ttTiJ 2t ~ ^' y -'- R ^ L 

^ L m,n 

where T^^" is the transmission matrix for an electron 
coming from the right. Electron energies are constrained 
by the energy conservation law, e+eL+ncu = e' +e R +mu>. 
Fermi functions fuR) describe the left (right) lead. We 
can write equation for current (^]) in a more compact 
form, by taking advantage of the time-reversal symme- 
try under which T R n ^{e' : e) -> T^(e, e') and the fact 
that the transmission matrix is nonzero when both elec- 
tron energies e and e' are within the band, 

Jtot = / deTT£$(e,e')[P(n)f L (e) (1 - f R (e')) 



e 

itTi 



n,m 



P(m)f R (e')(l-f L (e))} 



(9) 



We can also compute the elastic contribution to the total 
current in Eq. (g) by imposing the constraint of elastic 
tunneling n = m 



/It 
deVr!«fi(^ + A/i)P(n) 
-2* „ 



[fL(e)-f R (e + An)], 



(10) 



where A/x = — e R . The total energy of the electron 
E = e + EL = e'+e R is in this case conserved. In the limit 
of zero temperature Eq. ( JlOj ) gives the correct result for 
elastic tunneling (see e.g. Q). 

Our derivation of equation (^ is based on the one- 
electron approximation which leads to neglecting many- 
body effects as are a) the Coulomb repulsion, b) phonon 
mediated electron-electron interaction. Closely con- 
nected to the latter is our assumption, that phonon dis- 
tribution P(n) is independent of electron current through 
the tunneling region. Validity of the equation ^ is there- 
fore limited to cases, when the current is small and the 
tunneling region is strongly coupled to the external heat 
bath. This can be achieved even at large voltage bias 
when coupling to the leads is weak, i.e. to -C t. 



III. RESULTS FOR A SINGLE SITE 

In our calculation we have set the hopping within the 
lead to t — 1. The maximum number of allowed phonon 
quanta N p h was selected such that results fully converged 
for all chosen parameters of the system at different tem- 
peratures. 
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A. Transmission 

In Fig. (||) we show the total and elastic transmission 
(Eqs. (§|,0)) as a function of the incoming electron en- 
ergy e at different temperatures. At low temperatures 
we see a central peak positioned a,t e ^ e$ — \ 2 /u : a, one 
phonon side peak at e ~ eo + oj — \ 2 /u>, and a small 
two phonon peak at higher e. Positions of the peaks 
approximately correspond to polaron energies, given by 
e po i(n) = eo — \ 2 /oj+nu>; where, n represents the n— th ex- 
cited state of a polaron. By polaron we symbolize a state 
of an electron, coupled to phonon degrees of freedom lo- 
cated on the same site in the limit when to — > 0. It is im- 
portant to emphasize that all peaks have an elastic and 
an inelastic contribution (see also ||[l5|]). Elastic con- 
tributions correspond to tunneling through the ground 
or the excited state of a polaron without emitting (or 
absorbing) a phonon. Inelastic parts correspond to tun- 
neling through a given state with simultaneous phonon 
emission or absorption. At zero temperature only phonon 
emission processes are allowed since before tunneling the 
phonon state contains no phonons. 



in our numerically exact approach at small temperatures, 
and it changes at most 5% at large temperatures (see val- 
ues in the inset). We would like to stress, that sum-rules 
I j. (and I j see next subsection), defined in the work of 
Wingreen ||, do not represent sum-rules in a strict sense. 
They instead represent identities that are derived on the 
basis of two main approximations: dispersionless bands 
and a large bias limit used when calculating the current. 
Even though our calculations are not limited by these ap- 
proximations, we nevertheless chose to compare integrals 
of transmission with 1^, and in the next chapter integrals 
of current with Ij. As the temperature rises, peaks due 
to multi-phonon processes increase in strength. It almost 
seems as though finite temperature increases the effective 
coupling strength A However, an increased coupling 
strength would also shift the peaks. Peaks in Fig. (||) 
do not shift with increasing temperature. The increased 
strength of multi-phonon side-peaks is a consequence of 
increased number of excited phonons in the tunneling re- 
gion at higher temperatures. 

B. Current through a single site 
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We continue with a discussion of current. The relative 
positioning of bands and the choice of chemical potentials 
in the calculation is presented schematically in Fig. (|^). 
We have positioned the chemical potentials in the left 
and right lead in the middle of each band. This situa- 
tion corresponds to both leads being made from the same 
metal. The value of the chemical potential within each 
lead is constant and corresponds to asymptotic values in 
the leads. Bands are shifted symmetrically due to an ap- 
plied bias A/i i.e. ex, — A/x/2 and en = — A/i/2, thus 
ex, — £_r = A/i = eV; where, V represents the potential 
drop across the tunneling region. 



FIG. 2. Transmission probability as a function of the inci- 
dent electron energy calculated for different temperatures T. 
The heavy line is the total transmission Ttot, and the dashed 
(lower) line the elastic part T e i ast . The parameters of the 
Hamiltonian are: A = 0.25, q = e r — eo = 0.0 (no voltage 
drop across the dot), to = 0.2, ui — 0.5 and N p u = 17. Also 
presented are sum-rules where 7y = 2nt(to/t) 2 . 

At finite temperature other processes can take place. 
As an electron enters the tunneling region there may be 
one or more phonons excited in the system. The side 
peak at e ~ eo— o-> — X 2 /uj is due to a process when an elec- 
tron enters the tunneling region and absorbs a phonon. 
Such a process has been observed by Cai et al O]. As 
the temperature increases more inelastic channels open 
for electron tunneling, giving rise to an increased strength 
of side peaks. Interestingly enough, the sum-rule, valid 
at zero temperature and within the wide-band approxi- 
mation d|| It — I deTtot(e) = 2nt(to/t) 2 , remains valid 




FIG. 3. A schematic representation of the positioning 
of the bands. Some elastic and inelastic processes that con- 
tribute to the total current are shown. Shaded areas represent 
filled one-electron states. In elastic processes ({1} and {2}) 
the total electron energy is conserved e + el = e' + er, while 
in inelastic ones ({3} and {4}) it is not. For clarity, only zero- 
and one-phonon states on the central site are included. 

In Fig. (J4j) we present the total and elastic current vs. 
e (the on-site energy of the site j — 0) for four diffcr- 
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ent values of A/x. At a small A/x — 0.02 <C oj there is 
no inelastic contribution to the current (solid and dashed 
curves overlap). The main reason for this interesting ef- 
fect is that the chemical potential difference across the 
tunneling region is smaller than the minimum energy 
change u> necessary for the inelastic tunneling. The nec- 
essary condition for an electron, with an incoming en- 
ergy e, to contribute to the current at T — 0, besides 
having a finite tunneling rate, can be expressed with two 
inequalities: e + < /x^, and e' + £r + muj > [ir. At 
small Afi this condition leads to well defined peaks in the 
Jtot.eiasti^o) curves since only electrons with a total en- 
ergy in a narrow interval between the left and right chem- 
ical potential can give rise to current. Therefore only pro- 
cesses labeled by {1} and {2} in Fig. ((H) contribute to 
formation of the side (a) and the central peak (b) respec- 
tively. We should stress that phonons nevertheless still 
play an important role even at small voltage biases. The 
side-peak (a) at approximately eo = X 2 /lo — u> represents 
an elastic process where an electron elastically tunnels 
through the first excited state of a polaron whose energy 
is e po /(l) = — X 2 /uj + u>, while the main peak (b) corre- 
sponds to elastic tunneling through the polaron ground 
state e po i(0). 
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FIG. 4. The total current Jtot (solid line) and elastic cur- 
rent Jeiast (dashed line) vs. eo at T — for four different 
choices of A/x. All curves are in units of Ij (see definition be- 
low). Note also that the ordinates in figures A) and B) are dif- 
ferent than in C) and D). Current flows at different eo through 
different channels. Opening or closing of different channels 
(processes) as eo increases is reflected in peaks, shoulders or 
dips in the curves. Those anomalies (labeled by small let- 
ters) at different A^i correspond to the following processes (see 
Fig. (|)) A) a:{l}, b:{2}, B) a:{l}, b:{2}, C) a:{l}, b:{l,2,3}, 
b':{2}, c:{2,4}, D) a:{l}, b:{l,2,3}, c:{l,2,3,4} c':{2,4}. 

As A/x increases to A/x = u> = 0.5 the side peak de- 
velops into a shoulder and moves towards lower eo ~ 
lu 2 /X — ijj — A/i/2 while the central peak broadens. This 
is a consequence of separating the left and right chemical 



potentials and thus imposing less restrictive conditions 
on the tunneling electron energies. In particular, the 
broadening of the central peak is caused by simultane- 
ous tunneling through processes labeled {1} and {2} in 
Fig. (|^); also note the figure caption. Inelastic contribu- 
tion is still small. 

The main features appearing at larger bias, i.e A/i = 
0.8, are a) the emergence of the inelastic current (the 
solid and dashed lines do not overlap) and b) devel- 
opment of a new structure in the central peak region. 
The appearance of the inelastic current is caused by the 
opening of the inelastic channels at A/x > u>. The new 
structure in the central peak region is caused by open- 
ing and closing of elastic or inelastic tunneling channels 
as e increases [see caption to Fig. (|J)]. At even larger 
bias A/i — Aoj — 2.0 the current Jtot(^o) broadens even 
though the effects of different channel contributions are 
still visible. The inelastic current increases relatively to 
the elastic current. 

We have also computed integrals Ij = J deoJ(eo) 
that according to Wingreen et al. || should equal Ij = 
(e/irh)2irt(to/t) 2 An. Our findings are that even though 
the current does scale with A/x, nevertheless, our inte- 
grals deviate from sum-rules that exist in the wide band 
and large bias limit. The main reason is that the cur- 
rent flows through only selected channels allowed by the 
difference in the chemical potentials and not through all 
the channels as is the case in the limit of large bias and 
wide-bands. 

The effect of finite temperatures on the total current is 
presented in Fig. (||) . In contrast to the case of the total 
transmission we were unable to obtain additional peaks 
corresponding to phonon absorption processes. With in- 
creasing temperature the smearing by the Fermi func- 
tions overcomes development of those peaks. The effects 
of raising the temperature are more intense at smaller 
biases where peaks are narrower. Phonon side-peaks dis- 
appear around T ~ uj/2. 
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FIG. 5. The total current J tot vs. e at T = 0, 0.075, and 
0.2 for four different choices of Afi. All curves are in units of 
I j. All other parameters are the same as in Fig. (ph. 



IV. TUNNELING THROUGH 
AHARONOV-BOHM RING 

Next, we consider tunneling and consequently the cur- 
rent through an Aharonov-Bohm ring. The purpose of 
this section is to investigate the effect of dephasing by 
phonon modes on an electron as it tunnels through a 
ring. In detail, we will focus on the effect of dephasing 
on the transmission and current through the ring. 




FIG. 6. Schematic representation of the Hamiltonian de- 
scribing the Aharonov-Bohm ring. Dots surrounded by circles 
represents sites coupled to Einstein phonons. For simplicity 
we have chosen on-site energies on the ring to be constant 
ej — eo. Magnetic flux $ penetrates the center of the ring. 

The model is schematically presented in Fig. (|J). The 
ring consists of four sites, connected by hopping matrix 
element t' , that are coupled to two connecting leads with 
the hopping matrix element to- Each site of the ring is 
coupled to an Einstein phonon mode with frequency u>. 
There are four different phonon modes (one on each site) 
with identical frequency to. A magnetic flux is penetrat- 
ing the circle. This is reflected in an additional phase 
±<f> that the electron gains each time it hops from one 
site on the circle to another. The electron part of the 
Hamiltonian (Q), describing hopping within the ring, has 
to be modified in order to encompass the effect of a vector 
potential 

H el = e cfa ( e# 4 c i+i + h - c ) ' ( n ) 

3 3 

where the sum runs over the sites of the ring. All on-site 
energies within the ring were, for simplicity, set to eo. 
We solve this problem using the method described above 
(see also |lq]), and generalized to many phonon degrees 
of freedom. By increasing the maximum allowed number 
of phonons N p h on each site this problem can be solved 
with any desired accuracy. There are of course computer 
limitations. In practice, for moderate electron-phonon 



coupling and low temperatures, N p h = 3 is enough to 
obtain results with at least 1% accuracy. The number 
of many-body states increases as N st = N(N p f l + l) N = 
1024, which corresponds to the number of channels N c h = 
2N st , N is the number of sites in the ring. In a scattering 
problem where the energy of the incoming electron is 
known in advance, we have to solve a large (N st x N st ) 
sparse system of complex linear equations for each e. Due 
to computer limitations, we restricted calculations to zero 
temperature. 
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FIG. 7. Transmission probability as a function of the in- 
cident electron energy e at T — calculated at different 6. 
The heavy line is the total transmission T to t, and the dashed 
(lower) line the elastic part T e i as t- The parameters of the 
Hamiltonian are: A = 0.2, e; = e r — eo = 0.0, to = 0.3, 
u>o ~ 0.4 t' = 0.5 and N p h = 3. In the inset we show (insets 
are in the same units as the rest of the figures in Fig. (Q)) 
transmission T e i as t for the noninteracting case, i.e. A = 0. 

We have computed the total and elastic transmission 
vs. e through an Aharonov-Bohm ring at different values 
of total phase, 8 = <j>N . In Fig. (0) we show results for 
the total and elastic transmission through a ring. Posi- 
tions of the main peaks are approximately located close 
to energies corresponding to solutions of a tight bind- 
ing problem on a N = 4 site ring with a phase 4>, i.e. 
E(k) = -2t'cos(k + 4>) for k = 2ni/N. One and two 
phonon side peaks (shoulders) are also visible at approx- 
imately e rit k = E(k) + nu>. The height and the width 
of the main peaks change substantially when the posi- 
tion of the main peak coincides with the position of the 
phonon-side peak (see 6 = tt/4). 

The most interesting result is found at 8 = ir which 
corresponds to a phase, where in the case of noninter- 
acting tunneling due to negative interference, transmis- 
sion is exactly zero for any e. At finite electron phonon 
coupling we obtain, for 8 = tt, a finite total transmis- 
sion. However, the elastic part of transmission remains 
exactly (down to numerical accuracy) zero. This result is 
surprising considering that elastic tunneling processes are 
possible when an electron first emits and then reabsorbs a 
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phonon. A phase shift that the electron undergoes while 
emitting a phonon is thus exactly canceled by the phase 
shift after reabsorbing a phonon. Tunneling at 8 = n 
is therefore exclusively due to inelastic processes. This 
finding is in accordance with the well accepted fact that 
impurity (elastic) scattering can not destroy Aharonov- 
Bohm oscillations ]l6| ], since it does not cause phase de- 
coherence. 
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FIG. 8. The total current J to t (solid line) and elastic 
current J e i aa t (dashed line) calculated at A/i = 0.1 vs. e at 
T = for four different choices of phase 6. All curves are in 
units of I j. Note that the solid and the dashed lines overlap 
indicating Jtot = Jeiast- As in the previous section we set 
€l = A/i/2 and €r = — A/i/2. The rest of the parameters 
remain the same as in Fig. (Q). In the inset we show for 
comparison the elastic current J e iast for the noninteracting 
case. Insets are in the same units as the rest of the Fig. (^). 

Furthermore, we have calculated current, Eqs. 
and ([To|), vs. eo through the Aharonov-Bohm ring at 
small bias, A/i = 0.1, zero temperature and different 9. 
These results are presented in Fig. (||). As we have al- 
ready shown in the previous section, only elastic current 
can flow at small bias. Nevertheless, effects of electron- 
phonon interaction are clearly visible. There is no net 
current at 9 — tt since only the elastic part of transmis- 
sion can contribute to the elastic current. We predict 
that at small bias A/Lt <§; ui Aharonov-Bohm oscillations 
through the mesoscopic Aharonov-Bohm ring should not 
diminish significantly due to electron-phonon coupling, 
where coupling is limited to optical phonons. 

Last, we present results for current through the 
Aharonov-Bohm ring at large bias, A/i = 0.5 > io. Re- 
sults are presented in Fig. (|^). Two most important ef- 
fects of larger bias are broadening of the peaks and the 
appearance of the inelastic current {Jtot — Jeiast)- As a 
cosnequence of dephasing by inelastic scattering we ob- 
serve finite inelastic current at 8 = tt. 
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FIG. 9. The total current J to t (solid line) and elastic 
current Jeiast (dashed line) calculated at A/i = 0.5 vs. eo at 
T = for four different choices of phase 9. The rest is the 
same as in Fig. (^). 

V. CONCLUSIONS 

In summary, we have extended a numerically exact 
method for inelastic tunneling, to calculate trans- 
mission through a single site coupled to phonons at finite 
temperatures. We have further proposed an approximate 
formula for current based on exact results for the trans- 
mission matrix, taking into account filled Fermi seas left 
and right from the tunneling region and left and right 
flowing currents. Presented formalism gives correct result 
in the limit of elastic current. Our approach is entirely 
based on a one-electron approximation. By investigat- 
ing tunneling through an Aharonov-Bohm ring coupled 
to optic phonons, we have obtained numerically exact re- 
sults for transmission at zero temperature. We have also 
computed the current through the ring. 

We highlight some of the important findings of this 
work: 

Transmission through a single site vs. incident electron 
energy at finite temperatures displays additional peaks 
due to phonon absorption processes. The sum rule It, 
valid in the wide band limit and zero temperatures, re- 
mains approximately obeyed in our exact approach even 
at finite temperatures. Deviations are within 5% even at 
temperatures T > lj. 

Current through a single site mimics transmission 
curve at small bias A/i <c uj and T = 0. At voltage 
bias A/i < lu, and zero temperature only current due 
to elastic processes is possible, since inelastic processes 
are excluded due to filled final states. With increasing 
temperature features due to inelastic processes disappear 
around T ~ uj/2. Phonon absorption peaks do not ap- 
pear in the current at finite temperatures since they are 
smeared by Fermi functions. 

Transmission through the Aharonov-Bohm ring con- 
sists only of inelastic transmission at 8 — tt. The lack of 
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elastic transmission is considered as evidence that elas- 
tic processes, even though they are a part of inelastic 
electron-phonon scattering, do not change the phase of 
the tunneling electron. We have thus presented numerical 
proof that elastic scattering by phonons docs not affect 
Aharononv-Bohm oscillations. Dcphasing occurs only 
through inelastic processes. This is the reason why inter- 
ference effects, measured at small bias A/z < u> across the 
tunneling region, display strong Aharonov-Bohm oscilla- 
tions by changing the flux through the ring. Calculation 
of current at larger bias A/i > uu shows that only inelas- 
tic current flows through the Aharonov Bohm ring in the 
case when 6 — it. The voltage drop over the region has 
to be larger than the minimal energy difference necessary 
for inelastic scattering i.e. A/z > u>. 

It would be interesting to study the tunneling through 
the Aharonov-Bohm ring at finite temperatures. Such 
computation would demand larger Hilbert space and 
more computational time. Our method also allows the 
implementation of the nonlinear electron-phonon inter- 
action and linear and nonlinear phonon-phonon interac- 
tions. It would be fascinating to investigate, how phonon- 
phonon interactions, leading to internal phonon dynam- 
ics, affect dephasing of the tunneling electrons through 
an Aharonov-Bohm ring. 
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